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Abstract 

It is the purpose of the proposed research to develop a digital computer simulation 
model to study the effects of an internal fault in a large permanent magnet ac synchronous 
motor. Permanent magnet motors are being considered as an alternative for ships with 
electric propulsion systems. In an electric propulsion system a large motor will be directly 
connected to a propulsion shaft. A windmilling shaft will continue to turn the rotor of the 
propulsion motor after the motor has been disconnected from its electrical power supply 
source. 



Following an internal electrical fault in a propulsion motor, it is expected that the 
motor will be disconnected from its electrical supply source. With the ship operating at or 
near rated speed following a casualty to the propulsion plant, the ship will coast down to a 
stop or until the crew takes action to stop the ship. A windmilling permanent magnet 
motor will generate a large enough internal voltage to continue to support large fault 
currents. 

This research will focus on the fault transient and the motor behavior during the 
time that the propulsion shaft is windmilling. Shorting the motor terminals will be 
considered as a means of reducing the power input into the fault. 

Thesis supervisor: Dr. James L. Kirtley, Jr. 

Title: Associate Professor of Electrical Engineering 



2 



Acknowledgments 



I wish to thank my advisor. Professor James L. Kirtley for his support and 
guidance, without whom I would not have been able to complete this work. I would also 
like to thank CAPT A1 Brown and LCDR Jeff Reed for their intellectual support and 
Professor A. Douglas Carmichael, my thesis reader. 

I would like also to thank my parents and my wife’s parents for their support and 
encouragement. Finally, I want to thank my wife, Maribeth for her support and assistance 
throughout my graduate work. 

I wish to dedicate this thesis to my children, Barbara and Steffen for understanding 
and putting up with me during the past three years. 



3 



Table of Contents 



Abstract 

Acknowledgments 

Table of Contents 

Chapter 1. Introduction 

1 . 1 Ship Propulsion Systems 

1 .2 Permanent Magnet Motors 

1.3 Ship Model 

14 Faults in Permanent Magnet Ship Propulsion Motors 

1 . 5 Research Approach 

Chapter 2. Conceptual Design 

2.1 Air gap size 

2.2 Magnet dimensions 

2.3 Determination of terminal current and machine rating 

2.4 Machine Reactances 

2.5 Internal Voltage 

2.6 Stator sizing 

2.7 Stator leakage reactance 

2.8 Losses and machine efficiency 

2.9 Back iron sizing 

2.10 Weight Calculations 

Chapter 3. Dynamic Models 

3.1 Two axis transformation 

3.2 Per Unit Scaling 

3.3 Permanent Magnet Motor Model 

3.4 Network model 

3.5 Fault Model 

3.6 Motor Load Model 

Chapter 4. Simulation Model 

4. 1 Discussion of the simulation program 

4.2 Fault simulation 

4.3 Simulations 

4.4 Fault Parameters 



„ 2 

..3 

...4 

..6 

..8 

.9 

11 

12 

13 

15 

19 

19 

21 

22 

23 

25 

26 

27 

28 

29 

30 

30 

32 

39 

39 

41 

44 

44 

47 

50 

51 



4 



Chapter 5. Conclusions and Results 53 

5.1 Simulation Results 53 

5.1.1 Simulation with motor terminals open 53 

5.1.2 Simulation with shorting of the motor terminals 59 

5.3 Suggestions for Future Research 65 

5.4 Conclusions 65 

References 67 

Appendix A. Motor Design Spreadsheet 69 

Appendix B. Node and Network Pre-Processor 71 

Appendix C. Load Flow Program 76 

Appendix C- 1 . Load flow calculation program 78 

Appendix C-2. Node Voltage Calculation 82 

Appendix C-3. Line Admittance Calculation 84 

Appendix D. Line Simulation Input File 85 

Appendix E. Synchronous Machine / Network Simulator 92 

Appendix E- 1 Motor Object 109 

Appendix E-2. Network Program 124 



5 



Chapter 1. Introduction 



The use of electric ship propulsion offers significant advantages in ship design, 
construction and operation. Placing electric propulsion motors as far back in the ship as 
possible serves to reduce long shaft lengths. Propulsion shafts often go through many 
compartments, creating design complications with shaft line component alignment and 
compartment arrangement. In addition, electric propulsion can reduce the number and 
type of prime movers in the ship. Propulsion power and shipboard electrical power can be 
derived from common prime movers. 

In electric propulsion ships, prime movers can be located anywhere in the. ship. 
This flexibility can improve survivability and make maintenance and shipboard 
arrangements easier [1]. 

Using electric drive can provide the ship with greater operational flexibility. 
Electric drive ships can operate the prime movers at their optimal speed and most efficient 
speed. By operating the prime movers at their most efficient speed, the ship’s fuel 
consumption can be lowered. This could decrease the frequency of refueling and increase 
the endurance of the ship. 

Electric power propulsion has been used in ships, commercial and military, for 
over 50 years. However, some drawbacks have been higher costs and greater weight than 
mechanical drive systems. Advanced systems using permanent magnets have the potential 
of being lighter than similar conventional systems. 
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The US Navy is developing modern electric drive propulsion systems for its ships. 
In some of the conceptual designs for these electric propulsion systems, the use of 
permanent magnet propulsion motors and ship service generators has been considered. 

An issue that needs to be investigated is the behavior of such machines during and 
following an electrical fault inside the motor or generator, where the machine supply 
breakers might not be able to prevent damage. Of concern are internal arcing faults that 
can result in very large currents and high localized temperatures [2], These large currents 
and high temperatures can cause extensive damage to insulation and conductors and more 
importantly to the permanent magnet themselves. 

The arcing process itself can cause significant damage to the machine. Arcs can 
compromise the electrical, as well as, the mechanical integrity of the machine. Localized 
damage caused by the electric arc, such as pits or hardened spots, can serve as crack 
initiators from which cracks that can be induced by machine vibrations can initiate and 
propagate [3]. 

After a permanent magnet generator or motor is disconnected from the rest of the 
system it will continue to generate an internal voltage, E a f, until the field comes to a stop. 

In wound field machines the field circuit breaker can be tripped simultaneously with the 
main circuit breakers thus essentially eliminating the field. Figure 1 shows a simple 
schematic of an ac synchronous machine. 
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Figure 1-1 Simplified Motor 
Model 

The internal voltage generated in the machine is proportional to the rotational 
speed of the field. In big machines, such as those used for propulsion or power 
generation, these internal voltages can be very large. The generated internal voltage can 
be large enough to continue the arcing process and cause further damage while the 
machine is coasting down. 

The purpose of this research is to develop some tools that can be used to 
investigate what happens to a large permanent magnet motor following an internal fault 
while the motor coasts down. 

1.1 Ship Propulsion Systems 

Current naval propulsion systems consist of multiple diesel engines or gas turbines 
coupled to one or more propulsion shafts through a set of reduction gears and clutches. 
Using two engines per shaft provides redundancy and continuity of propulsion. In the 
case of damage or maintenance to one engine while the ship is underway, propulsion can 
still be maintained on that shaft. 
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Although mechanical systems are very simple and reliable, they have some major 
disadvantages. Mechanical drive systems require separate prime movers for electrical 
power generation. Alignment between the shaft and the prime mover requires that prime 
movers be located as low in the ship as possible. This requirement results in large 
amounts of “lost volume” inside the ship and the superstructure for intake and exhaust 
trunks. Some of this volume can be recovered in ships with electric drive. 

Relatively light prime movers, such as gas turbines, can be located higher in the 
ship as electric drive does not require alignment between the propulsion motor and the 
prime mover. Conceivably, prime movers could be located in the superstructure. 

Locating the prime movers higher in the ship can minimize the arrangeable volume lost to 
long exhaust and intake trunks. Other advantages and drawbacks of these systems are 
discussed in [4]. 

Another possibility that electric propulsion offers is the capability to move the 
propulsion motors outside the hull of the ship. 

1.2 Permanent Magnet Motors 

The machine used in this research is a permanent magnet ac synchronous motor. 
This is essentially an ac synchronous motor in which the field windings have been replaced 
by permanent magnets. The analysis used for this machine is that used for wound field ac 
synchronous machines assuming that the machine is excited by a field current of constant 
value. 

Permanent magnet machines have a number of advantages over wound field 
machines. One of their most significant advantages is the elimination of the field windings 
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and the power losses generated in these windings. In addition, since no electrical 
connections are required to supply the field, this eliminates the need for slip rings and 
brushes. Therefore, eliminating the field windings can also result in smaller machines than 
wound field machines. This is especially important in marine applications that are volume 
limited where space is critical. 

Some of the limitations of permanent magnet machines come from the permanent 
magnets themselves. Excessive currents in the motor windings or excessive heat could 
result in demagnetization of the magnets. Other limitations and advantages of permanent 
magnet machines are discussed in detail in [5] and [6], 

The control aspects and the electrical power requirements of permanent magnet ac 
machines are not addressed in this research. Such aspects of permanent magnet and other 
electrical machines are addressed in references such as [7] and [8]. 

For the purposes of this research the permanent magnet motor will be assumed to 
have been operating in steady state at the time the fault is initiated. After the fault is 
initiated the motor will be disconnected from the rest of the network and allowed to coast 
down. This research will only address large, low speed motors such as those that would 
be used for ship propulsion. In this configuration the propulsion shaft will be directly 
connected to the rotor of the machine. 

The permanent magnet motor used in this research is a notional machine. A 
notional machine was used as this type of propulsion machinery is not yet in use in naval 
ships. Design and sizing calculations of this machine are discussed in chapter two. 
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One of the advantages of electric drive propulsion is the ability to directly drive a 
propulsion shaft without the need for reduction gears. It is desirable that ship propulsion 
motors operate at low speeds, so these machines will have a large number of poles. As 
the number of poles increases the physical size of the machine will increase. During the 
notional design phase an effort was made to keep the size and weight of the motor 
comparable with current propulsion machinery. For this research a flux concentrating 
permanent magnet machine was selected. The geometry of this machine is shown and 
discussed in more detail in chapter two. 

1.3 Ship Model 

Following the fault to the motor, after the main supply breakers to the motor are 
open, the ship is assumed to coast down to some fraction of its initial speed. During this 
coast down period and during the initial phases of the casualty it will be assumed that no 
action is taken by the ship’s crew to slow down or stop the shaft. The basis for this 
assumption is that during the initial phase of the casualty, the first indication available to 
the operators will be the tripping of the main supply breakers. Once the main supply 
breakers trip, there will exist an inherent time delay before action is taken to stop the 
affected shaft and take corrective actions. This delay is due to the time that it will take for 
ship’s personnel to evaluate, recognize and act to combat the casualty. 

Automatic protection systems, that operate when the breakers trip, could be used 
to minimize the time delay in stopping the shaft. A problem with these systems is that they 
could initiate protective action during spurious breaker trips. During these trips protection 
is not necessary and could be detrimental to the operation of the ship. 
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The mechanical energy supplied to the motor while the ship is coasting down is 
proportional to the speed of the shaft squared. It can be shown that for a given propeller 
the rotational speed of a windmilling shaft is proportional to the speed of the ship. As the 
ship coasts down the mechanical energy supplied to the motor will decrease as the ship 
slows down. 

1.4 Faults in Permanent Magnet Ship Propulsion Motors 

The focus of this research is investigating the performance of large permanent 
magnet motors such as those that would be used for ship propulsion. It can be expected 
that part of the electrical protection of these large motors will be main supply breakers. 
These breakers will disconnect the motor from its electrical power supply upon detection 
of a fault in the motor. 

Internal faults are of special interest due to the large amounts of energy that can be 
supplied to these motors by a free spinning shaft after the electrical supply breakers are 
tripped. Since the field of the machine is supplied by the permanent magnets, the motor 
with the free spinning shaft will behave like a generator supplying power to the fault. 
Referring to figure 1-1, the excitation voltage, E a f, is proportional to the flux produced by 

the field, <t>f, and the frequency of the field excitation, co [5] 

E af =Ko4 f (1-1) 

In permanent magnet machines a constant field flux, 4>f, is supplied by the 

permanent magnets. Therefore to stop this generator action, the shaft must be stopped. 
Stopping the shaft could be accomplished by using a mechanical shaft brake, changing the 
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pitch of the propeller in ships with controllable pitch propellers or stopping or slowing 
down the ship. Ships with more than one propulsion shaft can use the unaffected shaft to 
slow down or stop the ship. 

Other ways to slow down or stop permanent magnet machines are discussed in [9]. 
These methods are for unfaulted machines whose shafts are not being driven by an 
external source such as a ship’s windmilling propeller. The windmilling propeller and shaft 
will rotate at a speed that is proportional to the ship’s speed. 

1.5 Research Approach 

This research studied the dynamic behavior of permanent magnet ac synchronous 
motors following an internal fault. The motor studied was assumed to be directly 
connected to a ship’s propulsion shaft. Following the fault, the motor was disconnected 
from its electrical supply source and allowed to windmill as the ship coasted down to a 
fraction of its initial speed. 

Different from wound field machines, permanent magnet machines have a constant 
field flux supplied by the magnets. As the shaft windmills this constant field flux will 
generate an internal voltage in the machine. A sufficiently large internal voltage can 
continue to support an internal fault in the machine as the shaft windmills. 

To accomplish the goals of this research the following tasks were identified and 
performed: 

1. A notional permanent magnet motor having the performance requirements of a 
ship’s propulsion motor was derived. To do this task, a motor design spreadsheet was 
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developed. By specifying the principal motor requirements the motor parameters were 
estimated. 

2. A dynamic model of the permanent magnet motor, which incorporated the 
parameters of the notional design was derived. 

3. Models for the internal fault and of the ship were derived. 

4. The models of the motor, fault and ship were incorporated into a dynamic 
simulation program. 

5. The simulation was run and the results evaluated. 
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Chapter 2. Conceptual Design 



This chapter describes the procedure used to design the notional permanent 
magnet motor used in this research. The motor designed incorporates desired attributes of 
a motor for naval propulsion. The procedure described in this chapter is implemented into 
an Excel [10] spreadsheet. Appendix A. For the motor design, a set of requirements was 
established and some initial assumptions were made as to the geometry and operating 
parameters of this motor. The initial assumptions and performance requirements are 
summarized in Table 1. 



Table 1. Motor Specifications 



Number of Phases 


3 


Frequency (Hz) 


60 


Rotor speed (rpm) 


200 


Rated power (Hp) 


40000 


Operating voltage (V) 


1000 


Power factor 


0.8 


Winding factor, k w 


0.9 


L/D 


0.22 


Tooth fraction, A. p 


0.5 



A line frequency of 60 Hz was selected for the motor design since this is the 
frequency commonly used in U. S. Navy ships electric service applications. 
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The type of machine that was used in this research is a flux concentrating 
permanent magnet machine. In this type of machine the magnets are oriented so that their 




Figure 2-1. Motor cross section 



Using the geometry described in figure 2-1 and the requirements specified in table 
1, some of the initial sizing considerations were started. For the required rotor speed, n, 
and electrical frequency, f, the number of pole pairs, p, in the machine was determined 
using p = 60f/n. A rotor speed of 200 rpm was selected for this design because this speed 
is consistent with rated shaft speeds for naval ships. 

Since this machine will be used for ship propulsion, one of its desired attributes is 
that it is comparable is size and weight to conventional propulsion machinery. By 
calculating the radius and the length of the machine, an initial determination of the 
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feasibility of the motor can be obtained. Weight calculations are performed once all the 
machine components are sized. 

Once the rated speed of the rotor was established, the radius of the rotor, R, was 
determined from the mechanical power output required and the average gap shear 
stress, (x) . To calculate the radius of the rotor the following two equations were used 

Power = (T orque)© m (2- 1 ) 

Torque = (x)(2tcRL)R (2-2) 

The aspect ratio indicated in table 1 was calculated using the following relation 
suggested by Levi in [10] 



_L 

D 




(2-3) 



The aspect ratio estimated using equation (2-3) minimizes the winding resistance for a 
given electromotive force. 

Using the aspect ratio, L/D, calculated from (2-3) and combining equations (2-1) 
and (2-2), the radius of the rotor is 



R 



Power 



UmOmMfL/D) 



(2-4) 



Once the machine radius was calculated, the active length of the machine was determined 
as L = 2R(L/D). 



A value for the average gap shear was estimated in order to calculate the size of 



the machine. The average gap shear stress can be estimated by (x) ~ -^=B,K , where Bj 



is the peak value of the fundamental magnetic flux density and K is the root mean square 
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(rms) value of the effective armature reaction current density [6]. Gap shear has units of 
pressure, pascals, (Pa). 



For this machine design the value selected for the average gap shear was based on 
machine sizing considerations and cooling requirements. By maximizing the gap shear, the 
machine radius can be reduced, equation (2-4); however, high values of shear will require 
additional cooling requirements. 

The value of shear stress selected is that of an air-cooled machine. This type of 
machine was selected in an effort to maintain simplicity. For this type of machine, and for 
the specified rating, a shear value of approximately 50,000 Pa was obtained [6]. 

Figure 2-2 shows the variation of rotor radius as a function of gap shear. From 

equation (2-4) it can be noted that R x . Higher shear values will result in smaller 

machines. However, these smaller machines will require additional cooling provisions, 
adding to the complexity and most likely to the overall size of the machine. 



Fractional reduction in Rotor Size as 
a Function of average gap shear 

1 * 

0.8 
0.6 
0 4 
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Figure 2-2. Radius reduction as a function of gap shear 
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Considering the above reasons an average gap shear of 50,000 Pa was selected and 



used for the motor design. This shear value is at the high end for typical slow speed air- 
cooled machines [6]. For the specified shear and aspect ratio the rotor diameter was 
calculated to be approximately 4.2 m and the length of the machine, L, approximately one 
meter. 

2.1 Air gap size 

Mechanical considerations and performance requirements drive the length of the 
air gap, g. The minimum physical length, in meters, of the air gap can be estimated from 
the following relation [11] 



Using the calculated radius of the rotor, equation (2-4), and L/D, (2-3) a minimum 
gap length of 5 mm is obtained. This gap length was used for this research. 

2.2 Magnet dimensions 

Considering the geometry shown in figure 2- 1 , and assuming that the magnets 
occupy one half of the circumference of the rotor, the width of a magnet is given by 
7lR 0 

w m = — and from figure 2-1 w m = 2Rsin— L . Combining these two equations and 
P 2 

solving for the magnet angle, 9t 



g* 3.35 x 10 




D (2-5) 



9 t = 2 sin 1 — (2-6) 

2p 



Then the pole face angle, 9 m , is given by 
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(2-7) 




The remaining dimension of the magnets yet to be determined is the height, h m . 
This dimension is determined by considerations other than just geometry, such as the 
magnetic flux density at the air gap. In order to calculate the magnet height, the gap flux 
density needs to be determined. 

In flux concentrating machines, the flux density in the gap is greater than the flux 
density in the magnets. To calculate the fields in this machine a simple reluctance model 
presented in [12] was used. For the given geometry, the flux path in the rotor involves a 
magnet and one half of each of two adjacent pole pieces. For one pole piece the 
permeance of the air gap is given by 

R0 

P gap = fl 0 L 51 (2-8) 

g 

and the permeance of a magnet is 

P„ = (2-9) 

W„ 



The magnetic flux density in the magnets, B m , is calculated using a simplified 

magnetic circuit. This circuit consists of the magnet’s own permeance in series with one 
half of the permeance of each of two pole pieces. So that the flux density of the magnets 
is given by 



B = B 

m o 



P 



gap 



P gap + Pm 



and the flux density in the gap can be calculated 



( 2 - 10 ) 
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( 2 - 11 ) 



B = B — 

gap m r e 



Solving equation (2- 1 1 ) for the magnet height, the following relation is obtained 

(2-12) 



h _ l^g 



B 2 

m 



For a given radius and magnet angle, the ratio of gap flux to magnet flux will 
determine the magnet height. The gap flux density that can be achieved in the machine 
will be limited by stator teeth and back iron saturation. This limit will be checked later in 
the design. Once a magnet height is selected, the magnet spacing needs to be checked. 
For the given geometry the closest point between two magnets occurs at the interior 
corners. Using simple geometry the distance between two adjacent corners can be shown 
to be [12] 



s = 2(R - h m )sin — - w m cos— (2-13) 

2p 2p 

After checking that the magnet dimensions are compatible with the rotor size and 
the gap flux density has been selected, the surface current density, K, necessary for 

V2<x) 



operation of the machine at rated power can be calculated, K 



B, 



2.3 Determination of terminal current and machine rating 

For a machine with a small gap, it can be assumed that the magnetic flux is not a 
function of radial position. The space fundamental of this flux is then of the form 

B, = — sin^^-B (2-14) 

K l 
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With this value of flux, the internal voltage can be estimated using 




The terminal current into the machine, I t , is given by 




6N 



slots 11 slot 



(2-16) 



where N slots w slots = 27tRA p and K = J a h sIot k s . Once the internal voltage and the terminal 

current are known, the rating of the machine |P + jQ| = 3V t It can be determined if the ratio 
of terminal voltage to internal voltage, v = V,/E af , is known. A method to calculate v 

will be proposed later in the chapter; however, this method requires that the machine 
reactances be known. 

2.4 Machine Reactances 

The permanent magnets are in the main flux path of the armature. The presence of 
the magnets will make the machine salient since the direct and quadrature axes will be 
affected differently. Derivation of the direct and quadrature axes’ reactances is shown in 
several places; however, for consistency the notation used in [12] will be maintained. The 
direct and quadrature reactances, Xj and X q , are given by 




(2-17) 
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These formulas have a new unknown, the number of turns per phase in the stator, 
N s . The value of N s will be calculated once v and the rating of the machine are known. 



2.5 Internal Voltage 

To calculate the ratio of terminal voltage to internal voltage the phasor diagram 
shown in figure 2-3 is used. A two-axis representation of the machine is shown in figure 
2-3. The mathematical transformation used to derive the two-axis model is discussed in 
chapter three. 




Figure 2-3. Phasor diagram for Negatively Salient 

Motor 



From figure 2-3 the following set of phasor relations are derived 



E*=^V t 2 +(l t X q ) 2 -2V l IX q siny 



I d =I t sin(v|/+5) 



IX cosy 

5 = tan 1 

V + IX q COS V|7 



(2-19) 

(2-20) 

(2-21) 
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( 2 - 22 ) 



E* =E--(X,-X a )l j 

To solve this set of equations an iterative method is suggested in [12]. The end 
result of the proposed method is a value for the ratio of terminal to internal voltage, v. To 
implement the method, the set of equations (2-19) through (2-22) and the machine 
reactances, equations (2-17) and (2-18) are normalized with respect to the internal 
voltage, equation (2-13). After normalizing, per unit scaling, and some simplification the 
reactances are given by 



X ad = 






Pg 



B 



(2-23) 



Hs**, V2— 1 



aq Pg " P ' B, 



1 — sin - 

2 



(2-24) 



The set of phasor relations after normalizing and assuming operation at rated 



current ( i t = 1.0 p.u.) and power factor is given by 



e * = .Jv 2 +x aq -2x aq vsinv|/ 


(2-25) 


E X aq VC0S V 

o = tan - : 

v 4- x aq vsin V|/ 


(2-26) 


i d = sin(v|/ +5) 


(2-27) 


e af = e ’ ~( X aq - X ad)ld 


(2-28) 



Using the normalized reactances and the set of equations (2-25) through (2-28) the 
suggested method of [12] seeks a set of values that will solve the above set of equations in 
which the fixed point is e a f= 1.0. The iterative method is started by selecting an arbitrary 
value for v and solving equations (2-25) through (2-28) sequentially. For subsequent 
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iterations a new value of v new = v oid /y[e ^ is used and the iterative process is continued 

until e a f converges to 1.0. Once convergence is obtained the final value of v is used to 
determine the machine rating 

r\ 

3V t I t =-fr— R 2 IAKB,v (2-29) 

V2 p p 

If a value of terminal voltage is selected, table 1, then the machine terminal current 
can be calculated, equation (2-15). 



2.6 Stator sizing 

Once the terminal current is known the number of stator turns per phase is 
calculated, equation (2-16). The slot width and height, the number of slots and the size of 
the back iron need to be calculated to complete the initial sizing of the machine. If a tooth 
fraction of X p = 0.5 is used and assuming that R »g, then the slot width is calculated as 



27tR ... 

w, = (2-30) 

6Nk 



and the number of slots in the stator is given by 



N a „„ = 6N,k,X p 



(2-31) 



The slot height, h s , is determined principally by the insulation required by the 

armature windings. Commonly these requirements are such that the copper area in the 
slot is between 40% to 60% of the area slot [6]. For this design the slot copper fraction, 
kcm selected was 0.5. With these assumptions the slot height is calculated by 



h s)ot = 



K 



Ja k s k Cu 



(2-32) 
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Limits on the current density, J a , are established by the maximum temperature rise, 
9 W , in degrees Kelvin, above ambient temperature (40 C), allowed in the copper. The 

allowable temperature rise is based on the class of insulation used in the machine. For an 
air-cooled machine, the maximum current density based on thermal considerations can be 
estimated by 



J < .^wYcu 1 



K 



(2-33) 



where yc u is the conductivity of copper and h is the overall heat transfer coefficient. 
Equation (2-33) represents the energy balance between the heat generated in the copper 
and the heat removed by the cooling medium, air [11]. In equation (2-33) an overall heat 

watts 

transfer coefficient, h = 30- was used. This coefficient was derived for an air- 

K x meter" 

cooled machine and is based on empirical calculations discussed in [1 1]. 

Equations (2-1) through (2-32) are used in Appendix A to do the initial machine 
sizing. Based on the dimensions calculated by these equations and the established machine 
geometry, the weight of the machine was estimated. In addition, once the dimensions and 
geometry of the machine have been estimated, other machine parameters such as the 
machine efficiency and stator leakage reactance can be estimated. 



2.7 Stator leakage reactance 

The stator leakage reactance was calculated using the methods of reference [11]. 
Leakage reactance consists of (1) slot, (2) tooth top, (3) end winding, and (4) harmonic 
reactances. Tooth top leakage reactance is proportional to the gap length. In 
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synchronous machines with permanent magnet field the gap length is small and tooth top 
leakage reactance can be neglected [11]. 

For an initial estimate, the leakage reactance was assumed to consist of slot and 
harmonic leakage reactances. The slot leakage and harmonic reactances were calculated 
using the procedure outlined in [ 1 1 ]. The calculated reactance was increased by 1 0% to 
account for the effects of end winding reactance in the total leakage reactance. 

The leakage reactance was assumed to be the same for the direct and quadrature 
axes. This assumption is commonly made in the literature [11]. 

2.8 Losses and machine efficiency 

The power dissipated in the machine determines its efficiency. This dissipated 
power will determine the cooling and ventilation requirements for the machine to ensure 
that allowable temperature limits are not exceeded. The losses considered in Appendix A 
comprise: no-load losses in the iron, friction and windage losses, copper losses and load 
losses. 

The rotating magnetic field will result in heat losses in the iron due to eddy 
currents and hysteresis. These type of losses occur primarily in the stator teeth, the 
surface of the rotor poles and the structural parts of the machine exposed to alternating 
magnetic fields. These losses are proportional to the square of the peak value of B. 
Hysteresis losses are proportional to the frequency and eddy current losses are 
proportional to the square of the frequency. Core losses are estimated in Appendix A 
using the relations discussed in [1 1]. 
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Mechanical losses in the machine are the result of friction in the bearings and 
windage between the rotor and the stator. These losses were estimated using the relation 
presented in [1 1]. 

The copper losses were calculated by estimating the resistance of the armature 
windings and using the rated terminal current. To calculate the armature resistance, the 
mean length of conductor per phase in the armature was estimated as suggested in [1 1]. 
For a three phase machine the armature copper losses are 3R a I 2 , where R a is the armature 

winding resistance and I is the line current. 

Load losses are caused by the flux produced by the armature currents. These 
losses include eddy current losses in the support structures, pole surfaces and damper 
windings. The load losses were estimated as 1% of the armature copper losses [10]. 

The machine efficiency was estimated in Appendix A using 

Efficiency = 1 - *° SS ^- 
input 

where the losses are the sum of the individual losses described in the preceding 
paragraphs. The calculated machine efficiency was 98%, close to the efficiency predicted 
by [6]. 

2.9 Back iron sizing 

The thickness of the back iron, h c , must be sufficient to carry the machine flux. 
Assuming a sinusoidal air gap flux density, B gap , the back iron thickness is determined 
using 
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dS = 0 (2-34) 



S 



If the air gap flux is assumed to be constant, using the specified machine geometry 
and if R»g, equation (2-34) can be simplified to 



where B s represents the flux density in the back iron. 

The minimum back iron thickness is then calculated by using a back iron flux 
density close to saturation. A value of B s = 1.2 T, rms, was used for the machine design. 

2.10 Weight Calculations 

A machine weight was estimated by calculating a rotor weight and a stator weight 
for the given machine geometry. The material composition of the different components 
was assumed to be that of use in standard motor construction. The machine weight 
calculated does not include weight of foundations, the weight of a cooling system or the 
weight of an enclosure for the motor. 

The weight of the motor is expressed in long tons 1 (lton) for easier comparison 
with current or conceptual drives for naval ships. 




(2-35) 



One long ton = 1016 kilograms 
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Chapter 3. Dynamic Models 



To perform the fault simulations a model of the permanent magnet motor and of the 
electrical fault were derived. These two models were incorporated into a dynamic simulation 
model developed by Professor James L. Kirtley, MIT. This chapter describes the various 
models used in the simulation. Models for a permanent magnet ac synchronous motor, arc 
fault and interconnections between the motor and the power system are presented. 

3.1 Two axis transformation 

The permanent magnet motor model used in this research is based on the derivations 
presented in [1 1] through [14]. It assumes linear magnetics and sinusoidal stator winding 
distributions. Inductance and resistance values for the motor were calculated using the 
methods discussed in chapter two and Appendix A, 

Using the symmetry of cylindrical rotation, a coordinate transformation that accounts 
for the relative motion between the stator and the rotor is introduced. This transformation 
maps the stator winding variables to a reference frame that rotates with the rotor. In this 
reference frame, mutual inductances are independent of rotor position. Such transformation is 
commonly known as the Park’s transformation. A version of this transformation used in [5] is 
given by 
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(3-1) 



where 0 is some arbitrary angle. For a reference frame rotating with the rotor 0 = cot where © 
is the speed of the rotor. This transformation is orthogonal and power invariant, 
regardless of the choice of angle. The inverse transformation is given by 
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(3-2) 



If the Park’s transformation is applied to an arbitrary vector, F, representing a set of 
stator variables, such as currents or voltages, the stator variables are mapped to a fixed rotor 
reference frame. In this stationary reference frame the machine reactances are constant, 
independent of rotor position. 




= TF 



abc 



= T 'F 



dqO 



(3-3) 



The d and q subscripts are used to represent the direct and quadrature axes, 
respectively. The direct axis is aligned with the polar axis and the quadrature axis is aligned 
with the interpolar space. Zero-sequence variables are represented by the subscript 0 in 
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equations (3-3). The zero-sequence components represent components of armature currents 
that do not produce a net air gap flux [4]. 



3.2 Per Unit Scaling 



The models used in the simulation have been scaled to a per unit (p.u.) system. Scaling 
to a per unit system produces variables whose magnitudes are close to one. The base values 
selected for this scaling can be arbitrary. However, for this analysis it is convenient to use the 
motor’s rated power and voltage as the base values for the per unit scaling. Expressing the 
motor variables in a per unit system allows the comparison of this motor against other motors 
regardless of their ratings. 

For this motor the base quantities are defined by 

p b =|v b i b 

U-4) 

®o 

To convert the motor parameters to the per unit system, the ordinary variables were 
divided by the corresponding base quantities. A new quantity will be introduced into the model 
to change torque to the per unit system. This new constant is defined as the inertia constant, 

H; it represents the rotor’s kinetic energy at rated speed divided by the motor rated power. 

The inertia constant, H, has units of seconds and is expressed as 



H = 




(3-5) 



32 



In equation (3-5), J is the moment of inertia of the shaft and co m is the shaft’s rotational 
speed. In general the moment of inertia of a circular shaft can be calculated as J = J r pdV , 

V 

where p is the shaft’s material density [3]. It can be shown that for a circular shaft this is 
equivalent to J = Wk , where W is the mass of the shaft and k is its radius of gyration. For a 
solid circular shaft of radius, R, the radius of gyration is given by k = R/V2 . A value for H 
was calculated for this motor using the mass and rotor radius calculated in Appendix A. 

Since the rotor is directly connected to the propulsion shaft, the length of the shaft and 
the propeller needs to be included in the moment of inertia calculation. The weight of the 
propeller and of the propulsion shafts can be estimated from propeller and shaft weights from 
U.S. Navy ships or commercial ships. 

3.3 Permanent Magnet Motor Model 

The permanent magnet machine model is defined by a set of six coupled electro- 
mechanical equations. The electrical equations in this set are statements of Kirchoff and 
Faraday’s laws, which describe the voltage-current relations in the machine. The mechanical 
equations are statements of Newton’s Second Law of Motion. Coupling between the 
mechanical and electrical systems is through the dependence of electrical torque on flux 
linkages and through the dependence of flux linkages on torque angle, 5. 

The machine parameters used for the simulations are calculated in Appendix A using 
the methodology discussed in chapter two. A two-axis (d-q axis) representation of the 
synchronous machine is used. 
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In the d-q reference frame, the direct axis of the motor is represented by the following 
circuit based on the model presented in references [5] and [12]. 



Ra Xal 




Figure 3-1. Direct axis circuit representation 

where the constant current source represents the field generated by the permanent magnets. 

Similarly the quadrature axis of the motor can be represented by the following circuit. 



Ra Xal 




Xkq 



Rkq 



Figure 3-2. Quadrature axis circuit representation 

The following notation has been introduced in the above two models: 

X a d and X aq represent the magnetizing reactances, 

Xai represents the armature leakage reactance, which is assumed to be equal 
for the d and q axes, 

Ra is the armature resistance, 
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Xkd and Xkq are the damper winding reactances, and 

Rkd and Rk q are the damper winding resistances. 

Using the models shown in figures 3-2 and 3-3, the three phase ac synchronous machine is 
represented by the following set of equations 



V,=-^- + R a I d -V.fi> 



dt 



dX a 

V = — - + R I +1.& 

q dt 3 q d 

V — I D T 



dt 

dX 



L kd A kd 



v k ,=— ii+RJ 



dt 



kq kq 



(3-6) 



2 © m 

For permanent magnet excitation, the field is represented by the constant current source , If m . 



The flux equations for the machine are 



— L d I d + L ad I kd +L ad I fm 
^kd = k ad I d + L kd I kd + L ad I fm 

\ =L q I q+ L aq I kq 

^"kq k a qlq ^~L k qI k q 



where the inductances, Ld and L q , are defined by 

L d = L a , + L ad 
Lq = L a) + L aq 

The simulation code was written to be used with any ac synchronous machine, 
regardless of its rating. For this reason it is intended to be used with machines scaled in the per 
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unit system. To use the dynamic simulation code, the above set of equations were scaled in the 
per unit system. Dividing by the appropriate base quantities, the sets of equations (3-5) and (3- 
6) are represented in the per unit system by 
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(3-8) 



(3-9) 



(3-10) 



T e =(Vd i q- V fqid) 

The damper windings' voltages, v^d and v^q, were set to zero since these windings are 

normally shorted. Using a Thevenin equivalent circuit of the d-axis model, figure 3-1, the 
constant current source can be represented as a voltage source, E af = X ad I fm in series with the 

magnetizing impedance, X a d. 
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Figure 3-3. Thevenin equivalent of d-axis model 

With the model presented in figure 3-3, equation (3-7) can be rewritten as follows 

(3-11) 

By solving the above equations for the currents, the d-q model of the motor can then 
be incorporated into a simulation code that connects the motor to an external network. 

Solving equations (3-9) through (3-11) for the currents, ij and i q , and after some simplifications 
the following sixth order state-space representation of the motor was obtained 



dv d 

dt 

dv|' q 

dt 



+ - 3 - + cow +co 0 v d (3-12) 
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■®H'd+ W 0 V q ( 3 “ 13 ) 
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(3-14) 
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(3-15) 



(3-16) 
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(3-17) 



— — = — ^-(T - T ) 
dt 2H m e ' 



In equations (3-12) through (3-17) the following quantities are defined: 



x 

e" = — — vj/ kd = Voltage behind subtransient reactance, 

x kd 



x 

e d = -\|/ k = Voltage behind subtransient reactance, 

\ 



T d " = — — = D-axis open circuit subtransient time constant, 

®o r kd 



T" = — = Q-axis open circuit subtransient time constant, 



( 0 o r kd 



r ad = — — = Direct axis armature time constant, 
co r 



x' 

T — M 

*aq — 

®o r a 



= Quadrature axis time constant, 



x kd 



x d = x d = D-axis subtransient reactance, 



x" = x — = Q-axis subtransient reactance, 



L kq 



The stator currents, in the motor reference frame, have been defined by 



Vd-e" 



(3-18) 



V q +e? 

‘q = TT 



(3-19) 
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Equations (3-12) and (3-13) are the stator equations, which have time constants of the 
order of 1/co 0 = 0.0026 seconds. If the transients of interest are in the order of 0. 1 to several 

seconds then two stator equations, equations (3-12) and (3-13), will become algebraic 
equations provided that the following conditions hold 



Furthermore if co « co 0 then equations (3-12) and (3-13) simplify to 

Vd = v d (3 -12a) 

Vq - v q (3 - 1 3 a) 

With these two algebraic equations the machine model has been reduced to a fourth order 
model, equations (3-14) through (3-17). 

3.4 Network model 

The external network model used in the simulation code is based on the model derived 
in [1 5]. This model uses the machine currents to interface with the network. The definition of 
the network is accomplished by defining the nodes and branches of the network. This model is 
discussed in chapter four. 

3.5 Fault Model 

The fault model used for this research is an ac arc fault whose properties and 
characteristics are described in detail in [2]. Reference [2] describes the electric arc as a self- 
sustained electrical discharge having a low voltage drop and capable of supporting large 
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